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With the two flavor Nambu—Jona—Lasinio (NJL) model, we carried out a phenomenological study on the chiral 
phase structure, mesonic properties, and transport properties of momentum-space anisotropic quark matter. To 
calculate the transport coefficients we utilized the kinetic theory in the relaxation time approximation, where the 
momentum anisotropy is embedded in the estimation of both the distribution function and relaxation time. It 
was shown that an increase in the anisotropy parameter € may result in a catalysis of chiral symmetry breaking. 
The critical endpoint (CEP) is shifted to lower temperatures and larger quark chemical potentials as € increases, 
and the impact of momentum anisotropy on the CEP temperature is almost the same as that on the quark 
chemical potential of the CEP. The meson masses and the associated decay widths also exhibit a significant 
£ dependence. It was observed that the temperature behavior of the scaled shear viscosity 7/ T? and scaled 
electrical conductivity ce / T exhibited a similar dip structure, with the minima of both 7 / T? and Cel / T shifting 
toward higher temperatures with increasing €. Furthermore, we demonstrated that the Seebeck coefficient S 
decreases when the temperature rises and its sign is positive, indicating that the dominant carriers for converting 
the temperature gradient to the electric field are up-quarks. The Seebeck coefficient S is significantly enhanced 
with a large € for a temperature below the critical temperature. 
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I. INTRODUCTION 


The properties of strongly interacting matter described by 
the quantum chromodynamics (QCD) in extreme conditions 
of temperature T and density have aroused a plethora of ex- 
perimental studies in the last thirty years [1, 2]. The exper- 
iment studies performed at the Relativistic Heavy Ion Col- 
lider (RHIC) in BNL and the Large Hadron Collider (LHC) 
in CERN have revealed that a new deconfined state of mat- 
ter, the quark-gluon plasma (QGP), can be created at high 
temperature. Further, the non-central heavy-ion collisions 
produce the strongest magnetic fields and orbital angular 
momenta, which can induce a number of novel phenom- 
ena [3-5]. The lattice QCD calculation, which is a powerful 
gauge invariant approach to investigate the non-perturbative 
properties, has also confirmed that the phase transition is a 
smooth and continuous crossover for vanishing chemical po- 
tential [6, 7]. Owing to the so-called fermion sign prob- 
lem [8], lattice QCD simulation is limited to low finite den- 
sity [9, 10], even though several calculation techniques, such 
as the Taylor expansion [11, 12], analytic continuations from 
imaginary to real chemical potential [13, 14], and multi- 
parameter reweighting method [15], have been proposed to 
address this problem and improve the validity at high chem- 
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ical potential. More detailed reviews of lattice calculation 
can be found in Refs. [16, 17]. Alternatively, one also can 
rely on effective models, the Dyson—Schwinger equation ap- 
proach [18, 19] and the functional renormalization group ap- 
proach [20, 21], to study the chiral aspect of QCD for finite 
baryon chemical potential ug. Currently, there are various 
QCD inspired effective models, such as the Nambu—Jona— 
Lasinio (NJL) model [22, 23], Polyakov-loop enhanced NJL 
(PNJL) model [24, 25], Quark-Meson (QM) model [26, 27], 
and Polyakov QM (PQM) model [28, 29], which not only 
can successfully describe the spontaneous chiral symmetry 
breaking and restoration of QCD but also have been applied 
to explore the QCD phase structure and internal properties 
of the meson at arbitrary T' and ug. These model calcula- 
tions have predicted that at high chemical potential, the phase 
transition is a first-order phase transition, and with decreas- 
ing jp, the first-order phase transition has to end at a critical 
end point (CEP) and change into a crossover. At this CEP the 
phase transition is of second order. However, owing to vari- 
ous approximations adopted in the model calculations, there 
is no agreement on the existence and location of the CEP in 
the phase diagram. Furthermore, the rotation effects [30, 31], 
magnetic field effects [32-35], finite-volume effects [36—40], 
non-extensive effects [41, 42], external electric fields [43-— 
45], and chiral chemical potential effects [46-49] have also 
been considered in the effective models to provide a better 
insight into the phase transition of the realistic QCD plasma. 

Apart from the QCD phase structure information, the trans- 
port coefficients, characterizing the non-equilibrium dynami- 
cal evolution of QCD matter in heavy-ion collisions [50-52], 
have also attracted significant attention. The shear viscos- 
ity 7, which quantifies the rate of momentum transfer in a 


fluid with inhomogeneous flow velocity, has been success- 
fully used in the viscous relativistic hydrodynamic descrip- 
tion of the QGP bulk dynamics. The small shear viscos- 
ity to entropy density ratio ņ/s can be extracted from the 
elliptic flow data [53]. In the literature, there are various 
frameworks for estimating the 7) of strongly interacting mat- 
ter, e.g., the kinetic theory within the relaxation time approx- 
imation (RTA), QCD effective models [54-58], the quasi- 
particle model (QPM) [59, 60], and the lattice QCD simu- 
lation [61]. The electrical conductivity ce}, as the response 
of a medium to an applied electric field, has also attracted at- 
tention in high energy physics. The presence of ca not only 
can affect the duration and strength of magnetic fields [62], 
but also is directly proportional to the emissivity and produc- 
tion of soft photons [63, 64]. The thermal behavior of ce; has 
been estimated using different approaches, such as the micro- 
scopic transport models [65—67], lattice gauge theory simu- 
lation [68, 69], hadron resonance gas model [70, 71], quasi- 
particle models [72, 73], effective models [54, 74], string per- 
colation model [75], and holographic method [76]. Recently, 
studies of electrical conductivity in QGP in the presence of 
magnetic fields have also been performed [77-79]. Another 
less concerned but interesting coefficient is the Seebeck coef- 
ficient (also called thermopower). When a spatial gradient of 
temperature exists in a conducting medium, a corresponding 
electric field can arise and vice versa, which is the Seebeck ef- 
fect. When the electric current induced by an electric field can 
compensate with the current owing to the temperature gradi- 
ent, the thermal diffusion ends. Accordingly, the efficiency 
of converting the temperature gradient to an electric field in 
the open circuit condition is quantified by the Seebeck coeffi- 
cient S. In past years, the Seebeck effect has been extensively 
investigated in condensed matter physics. Very recently, the 
exploration has been extended to the hot QCD matter. For 
example, the Seebeck coefficient with and without magnetic 
fields has been studied in both the hadronic matter [80, 81] 
and QGP [82, 83]. In Ref. [84], the Seebeck coefficient has 
also been estimated based on the NJL model, where the spa- 
tial gradient of the quark chemical potential is considered in 
addition to the presence of a temperature gradient. 


In the initial stages of a HIC, the pressure gradient of the 
created fireball along the beam direction (denoted as longi- 
tudinal z direction) is significantly lower than that along the 
transverse direction. After the rapid expansion of the medium 
along the beam direction, the system becomes much colder 
in the beam direction than the transverse direction, which 
causes the QGP to possess a local momentum anisotropy, 
and this anisotropy can survive during the entire evolution 
of the medium [85]. In addition, the presence of a strong 
magnetic field can also induce a local anisotropy in the mo- 
mentum space. Inspired by the presence of momentum- 
space anisotropy in HICs, the primary objective of the present 
work is to study phenomenologically its effect on the chi- 
ral phase structure, mesonic properties, and transport coef- 
ficients in the SU(2) NJL model. To incorporate the mo- 
mentum anisotropy into numerical calculations, we follow 
the anisotropic distribution function parametrization method 
proposed by Romatschke and Strickland [86], f*"°°(p) > 


f'°(\/p? + €(p- n)?). Here, the isotropic momentum-space 
distribution function of particle (f°) is deformed by rescal- 
ing one preferred direction in the momentum space with the 
unit vector n and introducing the directional dependent pa- 
rameter €, which is used to quantify the degree of momentum- 
space anisotropy. This method has been extensively em- 
ployed to study phenomenologically the impacts of momen- 
tum anisotropy on various observables, such as photon pro- 
duction [87, 88], parton self-energy [89, 90], heavy-quark po- 
tential [91], and various transport coefficients [92-94]. The 
relativistic anisotropic hydrodynamics (aHydro) models can 
provide a more accurate description of non-equilibrium dy- 
namics compared to other hydrodynamical models [95]. As 
in most previous studies, the focus of the present work is 
on the weakly anisotropic medium (very close to the equi- 
librium) for which || < 1 and the distribution function can 
be expanded up to linear order in €. We found that even for 
small values of anisotropy, the effective quark mass and me- 
son masses change significantly compared to the equilibrium 
result. Unlike most momentum anisotropy studies of trans- 
port coefficients in QGP, in which the effect of momentum 
anisotropy is not considered in the different particle interac- 
tion channels, in the present work, we incorporated the mo- 
mentum anisotropy to the estimation of the relaxation time to 
better study the impact of € on transport properties of quark 
matter near the phase transition temperature region. 

This paper is organized as follows: Sect. II provides a brief 
review of the basic formalism of the 2-flavor NJL model. In 
Sect. III and Sect. IV, we present a brief derivation of the 
expressions associated with the constituent quark mass and 
meson mass spectrum in both an isotropic and anisotropic 
medium. Section V includes the detailed procedure for 
obtaining the formulae of momentum-anisotropy-dependent 
transport coefficients. In Sect. VI, we present the estimation 
of the relaxation time for (anti-)quarks. The numerical results 
for various observables are phenomenologically analyzed in 
Sect. VII. In Sect. VII, the present work is summarized with 
an outlook. The formulae for the squared matrix elements 
in different quark-(anti-)quark elastic scattering processes are 
presented in the Appendix. 


Il. THEORETICAL FRAME 


In this work, we start from the standard two-flavor NJL 
model, which is a purely fermionic theory owing to the ab- 
sence of all gluonic degrees of freedom. Accordingly, the 
lagrangian is given as [22] 


L = (id — Tro) + Gy) + (iyah 0) 


where y(7) stands for the quark (antiquark) field with two- 
flavors (u,d) Ne = 2 and three colors Ne = 3. Mọ de- 
notes the diagonal matrix of the current quark mass of up and 
down-quarks, and we take mo = m? = m4 to ensure isospin 
symmetry of the NJL lagrangian. G is the effective coupling 
strength of four-point fermion interaction in the scalar and 
pseudoscalar channels. 7 is the vector of the Pauli matrix in 
the isospin space. 


In the NJL model, under the mean field (or Hartree) ap- 
proximation [22, 23], the quark self-energy is momentum- 
independent and can be identified as the constituent quark 
mass mg, which acts as order parameter for characterizing 
the chiral phase transition. For an off-equilibrium system, 
the evolution of the space-time dependence of the constituent 
quark mass in the real time formalism can be obtained by 
solving the gap equation [96] 


Mq = Mo — 2GiTriS<(a, x), (2) 
where S<(x, y) = i(ab(y)v(a)) with x = (t,x) is the real 
time Green function in the coordinate space [97, 98], (...) 
denotes the thermal average, and the trace runs over spin, 
color, and flavor degrees of freedom. Transforming Eq. (2) to 
the phase space with the help of Winger transformation and 
introducing the quasiparticle approximation (see Ref. [99] for 
details), the gap equation can be further written as [96, 98— 
100] 


dp mg = 
ee Sp 
(27)? Ep falz, p) falz, p) ; 
(3) 
p° + m2(x) is the quasi-quark energy. As the 


Ma = Mmo +4NpNe J 


where Ep = 
NJL model is a non-renormalizable model owing to the point- 
like four fermion interaction in the lagrangian, an ultraviolet 
cutoff A is used to regularize the divergent integral. In the 
non-equilibrium case, the space-time evolution of the one- 
particle distribution function f(x, p) in Eq. (3) is described 
by the Boltzmann—Vlasov transport equation from the NJL 
model in the Hartree level [100-102]. By solving the Vlasov 
equation together with the gap equation concurrently, the con- 
stituent quark mass affecting the space-time dependence of 
f(a, p) can be determined self-consistently. 

To better understand the meson dynamics in HICs, it is 
useful to study the structure of the meson propagation in 
the medium. In the framework of the NJL model, mesons 
are bound states of quarks and antiquarks (collective modes), 
and the meson propagator can be constructed by calculating 
the quark—antiquark effective scattering amplitude within the 
random phase approximation (RPA) [23, 96, 98]. Following 
Refs. [96, 98], the explicit form for the pion (7) and the sigma 
meson (o) propagators in the RPA reads as 

2iG 

~ 1-2GT,,(z, ko, k) 
All information of a meson is contained in the irreducible 
one-loop pseudoscalar or scalar polarization function Iy, 
where the subscript M corresponds to pseudoscalar (7) or 
scalar (7) mesons. The space-time dependence of the po- 
larization function in a non-equilibrium system is given 
as [96, 98] 


Dm(z, ko, k) (4) 


Mar = Nee f TPT Ut - fale.) fale) 
M = INcLVe 4m Ep q\®,Pp q\t,Pp 
{1 (kô — k? — vir) | Ep + Ep-k 
2Ep-k kê — (Ep + Ep-k})? 
“ti al 
: (5) 
kô = (Ep 7 Ep-x)? 


Here, vz denotes the real value of the bound meson energy, 
and vm = 0 (2m,) for the m (o) meson. 


HI. CONSTITUENT QUARK AND MESON IN AN 
ISOTROPIC QUARK MATTER 


In an expanding system (e.g., the dynamical process of 
heavy-ion collisions), the space-time dependence in the 
phase-space distribution function is hidden in the space-time 
dependence of temperature and chemical potential. However, 
for a uniform temperature and chemical potential, i.e., for a 
system in global equilibrium, the distribution function is well 
defined and independent of space-time. Therefore, in the 
equilibrium (isotropic) state, to investigate the chiral phase 
transition and mesonic properties within the NJL model, one 
can employ the imaginary-time formalism. Actually, the re- 
sults in Ref. [96] have indicated that the real-time calculation 
of the closed time-path Green’s function reproduces exactly 
the finite temperature result of the NJL model obtained from 
the Matsubara’s temperature Green’s function in the thermo- 
dynamical equilibrium limit. In the following, we will briefly 
present the procedure for the derivation of the polarization 
function in the imaginary time formalism. In an equilib- 
rium system, m,, which is temperature- and quark chemical 
potential-dependent, can be directly calculated from the self- 
consistent gap equation in momentum space [22, 23]: 


dpm x 
Mg = Mo + 4GN:N, / Gaye E, (1 fe (P) F3(0)). 
(6) 


We can see that by taking the thermal equilibrium distribu- 
tion functions (Fermi—Dirac distributions), Eq. (3) has the 
same form as Eq. (6). The equilibrium distribution function 
of (anti-)quark fi@ is given by 


fog (P) = lexpl(Ep — Haq) 6) + 17°, (7) 


where 6 = 1/T is the inverse temperature of the system, and 
a uniform quark chemical potential py = Hua = —pMg.g 18 
assumed. It is noted that the ultraviolet divergence is not pre- 
sented in the integrand containing Fermi—Dirac distribution 
functions, so the momentum integral do not need to be regu- 
larized for finite temperatures. In the equilibrium, the meson 
propagator is given as [57] 


2iG 


Dm(ko, k) = 1 — 2GTI,,(ko, k)’ 


(8) 
where II,, at an arbitrary temperature and quark chemical 
potential is given by [103, 104] 


N: Ne = 
872 [((7mq mq)’ ko +k?) 


x Bo(ko, k, H, E; Mq) + 2A(u, T, Mq)]- (9) 


H (ko, k) = 


The minus (plus) sign refers to pseudoscalar (scalar) mesons. 
The function A, which relates to the one-fermion-line inte- 
gral, in the imaginary time formalism for finite temperatures 


and quark chemical potentials is given as [103] 


167? f d°p 1 

A= 3 2 expliwny) J EA 

(10) 
where wn = (2n + 1)7/ are the fermionic Matsubara fre- 
quencies and the sum of n runs over all positive and negative 
integer values. It is to be understood that the limit y — 0 is 
to be taken after the Matsubara summation. The function Bo 
in Eq. (9) relates to the two-fermion-line integral. At finite u 


and T’, Bo is defined as [103] 
Bo (im, k, H, T, Mg) 


16x? ; dp : 
= D,exrtinnt) es (27)? ((twn + u)? — E?) 


1 
x (liwn — iv + u)? — E2)? (11) 


where we have abbreviated E” = Ep-k = 4/ (p — k}? + m? 
and E = Ep for convenience, and after the Matsubara sum- 
mation on n is carried out, the complex frequencies 714 are 
analytically continued to their values on the real plane, i.e., 
im — ko — ie (e > 0) with ko being the zero component 
of the associated four-momentum. The full calculations of 
functions A and Bo at arbitrary values of temperature and 
chemical potential can be found in Ref. [105]. After evalu- 
ating the Matsubara summation by contour integration in the 
usual fashion [106], Eq. (10) is given as 


dp 
A=snt | aU + fe — 1), (12) 


where the Fermi—Dirac distribution function is introduced. 
Similar to the treatment of function A, after Matsubara sum- 
mation over n and the Matsubara frequencies iv; are analyti- 
cally continued to real values, Eq. (11) can be rewritten in the 
following form: 


Bo (ko, k, H, T; Mq) 


dp fy(p) + fe(p) -1 
= 161? je 4 SE 


X | 1 7 1 
(E+in)?—(E)? (E iv) — (2? 


Inserting Eqs. (12) and (13) to Eq. (9), we finally can obtain 
the expression for II), in the equilibrium state, which is for- 
mally the same as Eq. (5), except that the distribution func- 
tions are ideal Fermi—Dirac distribution functions rather than 
space-time dependent distribution functions. 


. (13) 


IV. CONSTITUENT QUARK AND MESON IN A WEAKLY 
ANISOTROPIC MEDIUM 


As mentioned in the introduction, the consideration of mo- 
mentum anisotropy induced by rapid expansion of the hot 
QCD medium for existing phenomenological applications is 
mostly achieved by parameterizing the associated isotropic 


distribution functions. To proceed with the numerical calcu- 
lation, a specific form of anisotropic (non-equilibrium) mo- 
mentum distribution function is required. In this work, we 
utilized the Romatschke—Strickland (RS) form [86] in which 
the system exhibits a spheroidal momentum anisotropy, and 
the anisotropic distribution was obtained from an arbitrary 
isotropic distribution function by removing the particle with 
a momentum component along the direction of anisotropy, 
fri (p) = fi8°(,/p2 + E(p- n)?), where n and £ are, re- 
spectively, the anisotropy direction and anisotropy parameter. 
As in previous studies [107], we shall restrict ourselves here 
to a plasma close to local thermal equilibrium, i.e., close to 
isotropy in the momentum space. Accordingly, the explicit 
form of anisotropic momentum distribution function in the 
local rest frame can be written as 


1 


exp[(/p? + (pn)? +m — ply)/T] +1 

(14) 
It is worth noting that for the anisotropic (non-equilibrium) 
matter, the T’ and u’ appearing in Eq. (14) lose the usual 
meaning of T and u in the equilibrium system and may 
have dimensionful scales related to the mean particle mo- 
mentum [108]. If we assume the system to be very close 
to the equilibrium (in the small anisotropy limit), then the 
parameters T” and p’ still could be taken to be the usual T 
and u, respectively, as performed in previous studies [92— 
94, 107]. The anisotropy parameter € is defined as é = 
(p3) /(2(pj))— 1, where p1 = |p—(p-n)-n| and pj = p:n 
are the momentum components of particles perpendicular and 
parallel to n, respectively. As the precise time evolution of € 
is still an open question, we assume that the anisotropy pa- 
rameter € in a local anisotropic system is constant and in- 
dependent of time, where —1 < € < 0 corresponds to a 
contraction of momentum distribution along the direction of 
anisotropy and € > 0 corresponds to a stretching of momen- 
tum distribution along the direction of anisotropy. In Eq. (14), 
the three-velocity of partons and anisotropy unit vector are se- 
lected as 


ie (p) — 


n = (sin x, 0, cos X), (15) 
p = p(sin 8 cos ¢, sin 0 sin ¢, cos 8), (16) 


where y is the angle between p and n, and p = |p| through- 
out the computations. With this choice, the spheroidally 
anisotropic term, ¿(n - p)? in Eq. (14) can be written as €(n- 
p)? = €p?(siny cos osin + cos x cos0)? = £c(0, ¢ġ, X). 
We further assume that n points along the beam (z) axis, i.e., 
n = (0,0, 1). It is essential to note that we shall restrict our- 
selves here to a plasma close to the equilibrium state and have 
small anisotropy around the equilibrium state. Therefore, in 
the weak anisotropy limit (|| < 1), one can expand Eq. (14) 
around the isotropic limit and retain only the leading order in 
€. Accordingly, the anisotropic momentum distribution func- 
tion in the local rest frame can be further written as [107] 


¿(n-p)? 
ZET 


‘mae = f° Pa- F’), (17) 


where the second term is the anisotropic correction to 
the equilibrium distribution, which is also related to the 
leading-order viscous correction to the equilibrium distribu- 
tion in viscous hydrodynamics. For a fluid expanding one- 
dimensionally along the direction n in the Navier-Stokes 
limit, the explicit relation is given as [109, 110] 


=F. (18) 


which indicates that the non-zero shear viscosity (finite mo- 
mentum relaxation rate) in an expanding system can also ex- 
plicitly lead to the presence of momentum-space anisotropy. 
At the RHIC energy with the critical temperature T, =~ 
160 MeV, r ~ 6 fm/c and 7/s = 1/42, we can obtain 
€ ~ 0.3. In principle, for the non-equilibrium dynam- 
ics of the chiral phase transition, a self-consistent numerical 
study must be performed by solving the Boltzmann—Vlasov 
transport equation together with gap equation in terms of 
the space-time dependent quark distribution as mentioned in 
Sect. II. However, because the non-equilibrium distribution 
function in the local rest frame of weakly anisotropic sys- 
tems has a specific form and the temperature and chemical 
potential appearing in Eq. (17) are still considered as free pa- 
rameters, the space-time evolution is not addressed. There- 
fore, in the non-equilibrium states possessing small momen- 
tum space anisotropy, just by solving the gap equation with 
anisotropic momentum distribution, i.e., Eq. (17), it is pos- 
sible for us to investigate phenomenologically the impact of 
momentum anisotropy on the temperature and quark potential 
dependence of the constituent quark mass. Accordingly, the 
gap equation, i.e., Eq. (2) or Eq. (6), is modified as 


co 2 2 
p-dp mq 0 0 pF, 
0 = |2N.MG 01 
| cf T? = Ia Sa +1 + SET 


+Mo9 — Mq. (19) 


Here, we have abbreviated F, = f9 (1 — f?) + f2(1— f$) for 
convenience. The momentum anisotropy is also embedded in 
the study of mesonic properties by substituting the anisotropic 
momentum distribution in the A function part of Eq. (5), thus 
obtaining 


2 2 
= p“dp| oo, eo _ 4 _ E&P“ Fp 
anal 7 [i+ ie OO 


In the € — 0 limit, the above equation reduces to Eq. (12). 
Similar to the treatment of function A, the weak momen- 
tum anisotropy effects can enter the Bo function (as given 
in Eq. (5) or Eq. (13)) by using the anisotropic momentum 
distribution. Without loss of generality, we selected the co- 
ordinate system in such a way that k is parallel to the z-axis, 
1.€., 

k = (0,0,k), |k|/=k. (21) 


We first discuss a simple case, i.e., k = 0, ko Æ 0. The 
computation of function Bo in a weakly anisotropic medium 


is trivial, viz, 


dp 
Bos =87° 
os =8T J Or E 


(f2"*°(p) + f2""°(p) — 1) 


1 1 
i (z — 2Eko + iesgn(ko) + kê + 2Eko — ar! ` 
(22) 


In the integrand of the above equation, there are two poles at 
E = Eo = tko/2if m < Eo. Applying the Cauchy formula 


= = Po + ind(x), (23) 


limeso 
T — E 


where P denotes the Cauchy principal value. The function 
Bo can finally be rewritten as 


2d f 
Bo =P | rl o(p) + fi) -1- SPF, 
ko q \*0 q \*0 ST e 


Mor m2@( E — mm? 24 
xy) (È) — mo? — m). (24) 
Here, zo = 4/ (42)2 — m2, and © is the step function to en- 


sure that the imaginary part appears only for kg/2 > mg. 
Next, in the case of k > 0, ko Æ 0, the expression for func- 
tion Bo is slightly complicated and can be written as 


°° pdp aani is 
TE (fa°(P) + JT (p) _ 1)) 


1 
1 
d 
à n a + (ke + 2koE — k?) /2pk — iesgn(ko) 
i ) 
x + (k2 — 2koE — k?)/2pk + iesgn(ko) 
= By + Bi, (25) 


+ 


iso 


with the abbreviation x = cos 0. The isotropic part B5° reads 


as 
iso 1 G Pp dp O ; g0 
Be =| f PAREN 1 


ali (ke — k? + 2pk)? — (2ko E}? 
(k2 + k2 — 2pk)2 — (2ko E)? 
+in(@(2pk — |k? + 2koE — k’ |) 


O(2pk — |k? — 2ko E eD]. (26) 


The anisotropic part Bọ” reads as 


(k2 — k? + 2ko E) 


M 1 [°° pd 2F 
panio Pp p Ep p ( 


k2 — kK? +2k)E\? 
(Be) me 


ke + 2koE — k? a 


kjo E 2ET pk 2pk kê + 2koE — k? — 2pk 
(kg — k? —2koB) | (kj — k? — 2koE F k2 — 2ko E — k? + 2pk 
pk i pk © | 2 — Qk B — k? — 2pk 
k2 + 2koE — k? k2 — 2ko BE — k? 
in( 0+ aa O(2pk — |k? + 2ko EB — k?|) — = sak O(2pk — |k? — 2ko E #))| . (27) 
p p 


When ko = 0, the imaginary part of Eq. (25) vanishes. There- 
fore, the real and imaginary parts of meson polarization func- 


tions for different cases in a weakly anisotropic medium are 
given by 


° dpp? 0 0 Ep? E? — vi,/4 
Relm (ko, 0) = NN: 1 0 F 28 
ulko) = NNP | PE] — Se) - 190) + Soh] (8) 
NN z2 
IT (ko, 0) = TENE fag — amg — R f1- Aleo) = Ao) + | O — 4d), 29) 
NN [© dpp? k? + v3, k — 2p o o 
Ilm = 1 l 1 y 
Remo) = Spt S a EP in| (1 — A) — 900) 
N.Nyz [°° Ep*dp k? +i, k? k — 2p 
1 | 
+f, er? 4 Apk |k-+ 2p (30) 


Similar to the treatment in Eq. (6), the vacuum divergence in 
Eqs. (28)-(30) also need to be regularized. When the effect of 
momentum anisotropy is turned off (€ = 0), Eqs. (28)-(30) 
are reduced to the results of Ref. [57] in thermal equilibrium. 
Once the propagators of mesons are given, their masses can 
be determined by the pole in Eq. (8) at zero three-momentum, 
i.e., 


1 — 2GRelIm (Mr o, 0) = 0. (81) 


The solution is a real value for My, < 2mg, and a meson 
is stable. However, for Mr o > 2m4, a meson dissociates to 
its constituents and becomes a resonant state. Accordingly, 
the polarization function is a complex function and II), has 
an imaginary part that is related to the decay width of the 
resonance as Tm = ImlIm (Mro, 0)/Mr,o- 


V. TRANSPORT COEFFICIENTS IN AN ANISOTROPIC 
QUARK MATTER 


In this section, we study the effects due to the local 
anisotropy of the plasma in the momentum space on the trans- 
port coefficients (shear viscosity, electrical conductivity, and 
Seebeck coefficient). The calculation is performed in the ki- 
netic theory that is widely used to describe the evolution of 
the non-equilibrium many-body system in the dilute limit. 
Assuming that the system has a slight deviation from the equi- 
librium, the relaxation time approximation (RTA) can be rea- 
sonably employed. The momentum anisotropy is encoded in 


the phase-space distribution function, which evolves accord- 
ing to the relativistic Boltzmann equation. We provide the 
following procedures for deriving the €-dependent transport 
coefficients. 


A. Shear viscosity 


The propagation of a single quasiparticle with temperature- 
and chemical potential-dependent mass in an anisotropic 
medium is described by the relativistic Boltzmann—Vlasov 
equation [101] 


1 
ph On + T mP falz, p) = Clfa(z,p)], (32) 


where 10t m? acts as the force term attributed to the residual 
mean field interaction. The right-hand side of Eq. (32) is the 
collision term. Considering that the system has a small depar- 
ture from the equilibrium due to an external perturbation, the 
collision term within the RTA can be given as 

Cif] ~ pP uulfalt, p) B f2(x,p)] = peupofa 


Ta a 


» (33) 


in which 7, denotes the relaxation time for particle species a 
and can quantify how fast the system reaches the equilibrium 
again. Near equilibrium, the distribution function can be ex- 
panded about a local equilibrium distribution function for the 


quarks as 

fa = fa (2: P) + Ofa, (34) 
with fa (2, p) = [exp( (ur (£)p” —Ha(x))8(x))+1]-*, where 
p= (Ex, p) is the particle four-momentum, and u” = 


7,(1, u) is the fluid four-velocity with y, = (1—u?)!/?. df, 
in Eq. (33) is the deviation of the distribution function from 
the local equilibrium due to the external disturbance, which 
up to first-order in the gradient expansion can be read as 


bf, =-—* 


pup 


pran fe +ma ae (O"T)O”) fI, (35) 


where the four-derivative can be decomposed into 0, = 
0/0" = u,D+V,, and D = ud, and V, = A#”d, 
denote the time derivative and spatial gradient operator in the 
local rest frame, respectively. g'” = diag(1,—1,—1,—1) is 
the metric tensor, and A#” = gt” — uu” is the projection 
operator orthogonal to u”. In the presence of weak momen- 
tum anisotropy, the associated covariant version of the weakly 
anisotropic function f$” for particle species a can be written 
as 


aniso C2 1 
©?) = ole EVE FI’ 
(36) 
~ fe (pt YY" oq ree (37) 


2T p’ Up 
where V” = (0,n) is defined as the anisotropy vector. Em- 
ploying Eq. (37) in Eq. (35), ô fa can be decomposed into two 
parts: 


Sfo = Ofa + bia (38) 


The first term on the right hand side is 


a(l — fe) pp” plu, Dur T V iti) 


T, 
ofie = a 
puy 


dm 
+p! (pM ty.) (ty. DB + Vab) + m prl, (39) 


Employing the equation of motion in ideal hydrodynamics 
and ideal thermodynamic relations, ô f'*° can be rewritten as 


i Ta ptp” 
fe — 0 1 0 y 
fie = pa- o, 
dm., > 

dT? x) 


1 6 
+| (Mun? i Auer" tko 


with 6 = Q,u®% and c? being the expansion rate of the 
fluid and squared sound velocity in the medium, respec- 
tively. The velocity stress tensor has the usual definition: 
aH” = LAKAY (V aug + Veta — 4A#”0). After tedious 


calculations, one can obtain the second term in Eq. (38): 


apre = VREE opo- nyape 
ees pp’ Blu, Duy +V pu) 
-p"(p"up)(uu DB + Y ub) — a pr] 
aAa) 
= VEE oe — nyag 
EEA 
C ETA 


(41) 


Allowing the system to be slightly out of equilibrium, the 
energy-momentum tensor T#” can be expanded as T#” = 
TE” +The Where T$” is the ideal perfect fluid form and TK 
is thie dissipative part of the energy-momentum tensor. In the 
hydrodynamical description of hot QCD matter, the dissipa- 
tive part of the energy-momentum tensor up to the first order 
in the gradient expansion has the following form [111]: 


Hv pv 
Tiiss = m 


— IAK”, (42) 


where 7” and II are the shear stress tensor and bulk viscous 
pressure, respectively. In present work, our focus is the shear 
viscosity component only. In the kinetic theory, the first-order 
shear stress tensor 7“” can be constructed in terms of the dis- 
tribution functions 


uv dp 1 uv p? p 
T - | h 7 wpe P of. (43) 


Here, the double projection operator is defined as Ay = 
AGA + APA — AY A, which can project any 
rank-2 Lorentz tensor onto its transverse (to u”) and traceless 
part. Inserting Eqs. (38)-(41) to Eq. (43) and comparing with 
the first-order Navier-Stokes equation 74” = 2ņo”” [112], 
in the rest frame of the thermal system with u” = (1,0) 
and p’u, = Ea, we finally obtain the expression of the €- 
dependent shear viscosity of particle species a, 


_ Eda [3 Tap? 0 0 T 
Na = 180T2 2 E3 fa fa) 2fa op ED 
i da dp Tap? 0 0 
" 30T T2 E2 fa fah (44) 


which is consistent with the result from Ref. [92]. For a sys- 
tem consisting of multiple particle species, the total shear vis- 
cosity is given as n = „Na. For SU(2) light quark matter, 
a = u, d,ū,d and the spin-color degeneracy factor reads ex- 
plicitly as da = 2N}. 


B. Electrical conductivity and Seebeck coefficient 


We also investigate the effect of momentum anisotropy on 
the electrical conductivity and thermoelectric coefficient. Un- 
der the RTA, the relativistic Boltzmann—Vlasov equation for 
the distribution function of single-quasiparticle of charge ea 
in the presence of an external electromagnetic field is given 
by 


1 p 


Ta 


(45) 


where F*”” is the electromagnetic field strength tensor. We 
only consider the presence of an external electric field, FP’? = 

—F® = € = (E,0,0). It is convenient to work in the local 
rest frame of plasma, and under the gedy state assumption 
(fa does not depend on time explicitly, 5#=0), Eq. (45) can 
be given by [84] 


Ofa _ Ofa 
(euE — V Ea) ý dp = ty! (46) 


Va Via F 


where we have used of the chain rule 22 8 p2 s2 e 
dp dp? ' dp dp’ 


is the electric charge of particle species a, and Va = OE, /Op 
is its velocity. To solve Eq. (46), we assume that the deviation 
of the distribution function in an anisotropic medium satisfies 
the following linear form: 


7 OEa oh ea ð aniso 
Of, = —Talea€ ax J- ap TaVa ` ae (47) 


The spatial gradient of the equilibrium isotropic distribution 
Ox f° in the presence of a medium-dependent quasiparticle 
mass can be expressed in the following linear form: 


0 0 0 Ea Ha 

Ox fa = fal fa) (oxi T ) ache) , (48) 
where Ha = tau denotes the chemical potential of particle 
species a and ta = +1(—1) for the quark (antiquark). Con- 
sidering u to be homogeneous in space and a temperature gra- 
dient only existing along the x-axis, and inserting Eq. (47) 
into Eq. (46), the perturbative term ôf, in an anisotropic 
medium can then be written as 


Of, = HaTa(eaEvx) — GaTaða( E va 
E(p-n)? ao 0 
E, T a(l fa)TaðrEa. (49) 


The expressions of H, and F, in the above equation are as 
follows, respectively: 


Ha = zf2(1— FNL + (8, 4,50) 


Ep’c(9, b, x) T 
JE T2 fall 18) (1-298 + 2) 450) 
o1 ép’c(0, b, x) 
Ga = alel fe) QE, TS (Ea — Ha) 
a- (1 2f? = i (51) 
Ea — Ha 


In the linear response theory, the general formula for the elec- 
tric current density Ja due to particle species a in response to 
an external electric field (£) and temperature gradient (V+T) 
is given by [113] 


Ja = Oela lE — SaVaT), (52) 


where Cela and Sa are the electrical conductivity and See- 
beck coefficient of the a-th particle, respectively. In terms of 
the distribution function, J, within the kinetic theory can be 


written as 
= cada [as sa Web fa (53) 


Finally, the expressions for Cea and Sa in the weakly 
anisotropic medium are respectively obtained as follows: 


e?daT, dp pt 0 E 
Oel a 6T [3 gzel #2) ( +$) 


Ee?daTa / dp p? pa f°) (1 _ 2f? + z) f 


36T? 1? E3 E, 
54 
and 
1 €adaTa dp pt i ò 
a= E, 7 {= 
a Oel,a | 6T2 I T2 Fz. H )fal fa) 
€ada€Ta dp pê 0 6 
Ea 2 J= 
3672 | 72 pa PeT Koal — fa) 
T a 
x {1-2 = a (55) 
( i (Ea =, =a) | Cela 


where a, is the thermoelectric conductivity due to particle 
species a. In the isotropic limit € — 0, Eqs. (54)-(55) 
reduce to the formulae in the equilibrium. In condensed 
physics, a semiconductor can exhibit either electron con- 
duction (negative thermopower) or hole conduction (positive 
thermopower). The total thermopower in a material with dif- 
ferent carrier types is given by the sum of these two contri- 
butions weighted by their respective electrical conductivity 
values [114, 115]. Inspired by this, the total Seebeck coeffi- 
cient in a medium composed of light quarks and antiquarks 
can be given as 

S= Da SaCel,a = A Qa — Q i (56) 
yx Cela D Oel,a Oel 


where the fractional electric charges of up and down (anti- 
)quarks are given explicitly by eu, = —ea = 2e/3 and eg = 
—ejz = —e/3, where the electric charge reads e = (47a,)!/? 
with the fine structure constant a, ~ 1/137. 


VI. COMPUTATION OF THE RELAXATION TIME 


To quantify the transport coefficients, one needs to spec- 
ify the relaxation time. In present work, the scattering pro- 
cesses of (anti)quarks through the exchange of mesons are 
encoded into the estimation of the relaxation time. The relax- 
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Fig. 1. (Color online) (a) Temperature dependence of the con- 


stituent quark mass mg at 4 = 0 MeV for different fixed anisotropy 
parameters. (b) Chiral susceptibility Xen at y = 0 MeV for different 
fixed anisotropy parameters. The broad dashed lines, dashed lines, 
and solid lines correspond to the results for € = —0.3, € = 0, and 
0.3, respectively. 


ation times of (anti)quarks are microscopically determined by 
the thermal-averaged elastic scattering cross section and par- 
ticle density. For light quarks, the relaxation time in the RTA 
can be written as [57] 


Tp (T, y) = Nqluasua + Cugsad + udud 
+Ng [Fud—ud + TETN (57) 


where the number density of (anti-)quarks in a weakly 


à : Š : : = da’ p fan 
anisotropic medium is given as nyg = df tmp Jaa) 


with d; = d; = 2N, denoting the degeneracy factor. The 
momenta of the colliding particles for the elastic scattering 
process a(p;) + b(p2) — c(p3) + d(p4) obey the rela- 
tion pı + P2 = p3 + p4 = 0, and we use the notation 


|pi| = |p2| = p for convenience. In the center-of-mass (cm) 
frame, the Mandelstam variables s, t, u are given by 


s= 4m? + 4p”, t = —2p?(1 — cos p), 
u = —2p°(1 + cos 6p), (58) 


where 0, is the scattering angle in the c.m. frame. The Man- 
delstam variables hold the relation u + s + t = 4m. Cab cd 
denotes the thermal-averaged elastic scattering cross section 
in the weakly anisotropic system, which can be written as 


co tmax da b d 1 
Baa = f ds f m sin? o | dx3 
So t ao} 
1 


min 


(1 — F2 (pem, Hy 23)) / dra 


=1 
x(1— fF (Pem, H, x4))L(s, u, £1, £2), (59) 


with (1 — fy) denoting the Pauli-blocking factor for the 
fermions due to the fact that some of the final states are 
da 


already occupied by other identical (anti-)quarks. % = 


Torse mzy M |? is the differential scattering cross section, 
with |/|? denoting the squared matrix element of a specific 
scattering process. The formulae of the squared matrix el- 
ements for various scattering processes are presented in the 


Appendix. The integration limits of t are tmax = 0 and 
tmin = —4P2n = —(8 — 4m2), with Pem = (4/8 — 4m?) /2 
denoting the momentum in the c.m. frame. The kinematic 
boundary of s reads sg = 4m. The scattering weighting fac- 
—4t(s+t—4m?) 
(s—4m2)? 
tering processes with a small initial angle because the large 
angle scattering dominates in the momentum transport pro- 
cess [116]. In the c.m. frame, the leading-order anisotropic 
distribution function can be rewritten as 


tor sin? 6, = is introduced to exclude the scat- 


FP (Pem, H, x) = f? (Pem, Lt) 


2 2 
— Fink (Pam p)(L = f (Pem )), (60) 


—m? m? 
oe = \/s/2. In Eq. (59), £ denotes the 
probability of finding a quark-(anti)quark pair with the center 
of mass energy ,/s in the anisotropic medium, and it is given 
by 


where Eon = 


aniso 


L =C s(s _ 4m2) Fa (Pem, H, xı) b (Pem, H, x2) 
x Vrei (8), (61) 


[s—4m2 . f . 
where U;e1(S) = SMa is the relative velocity between 


two incoming particles ‘in the c.m. frame; C is the normal- 
ization constant, which is determined from the requirement 
that i ee ds J} dz J} deol = 1, where x; = cos, and 
z2 = cos 02, with 6; being the angle between p; and n and 
zı = —X3, L2 = —x4. Applying the above formula of re- 
laxation time to Eqs. (44), (54), and (55), we can calculate 
the transport coefficients in the QCD medium and study their 
sensitivity to the momentum anisotropy. 
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Fig. 2. (Color online) (a) Double mass of the constituent quarks 
2mq (gray lines), 7 meson mass (blue lines), and o meson mass 
(red lines) as a function of temperature at u = 0 GeV for different 
anisotropy parameters €. (b) Temperature dependences of m (blue 
lines) and o (red lines) meson decay widths for different €. The 
long dashed lines, dashed lines, and solid lines represent the results 
for € = —0.3, € = 0, and 0.3, respectively, with the corresponding 
Mott temperatures approximately given by 187 MeV, 196 MeV, and 
206 MeV. 


VII. RESULTS AND DISCUSSION 


Throughout this work, the following parameter set is used: 
mo MO, u ™0o,d 5.6 MeV, GA? = 2.44 and A = 
587.9 MeV. These values are taken from Ref. [117], where 
these parameters are determined by fitting quantities in the 
vacuum (T = u = 0 MeV). At T = 0, the chiral symmetry is 
spontaneously broken and one obtains the current pion mass 
Mo,» = 135 MeV, pion decay constant f+ = 92.4 MeV, and 
quark condensate — (77)!/3 = 241 MeV. 

In the NJL model, the constituent quark mass is a good 
indicator and an order parameter for analyzing the dy- 
namical feature of chiral phase transition. In the asymp- 
totic expansion-driven momentum anisotropic system, the 
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anisotropy parameter € is always positive owing to the rapid 
expansion along the beam direction. However, in the presence 
of a strong magnetic €, it becomes negative because of the 
reduction in transverse momentum due to Landau quantiza- 
tion. As we restrict the analysis to only a weakly anisotropic 
medium, the anisotropy parameter we address here is artifi- 
cially taken as € = —0.3, 0.0, 0.3 to investigate phenomeno- 
logically the effect of € on various quantities. In Fig. 1 (a), 
we show the thermal behavior of the light constituent quark 
mass mq for vanishing quark chemical potential at different £. 
For low temperature, mą remains approximately constant at 
(mq ~ 400 MeV), and then, with increasing temperature mg, 
it continuously drops to near zero. The transition to small 
mass occurs at higher temperature for a higher value of £. 
These phenomena imply that at zero chemical potential, the 
restoration of the chiral symmetry (the chiral symmetry is not 
strictly restored because the current quark mass is nonzero) in 
an (an-)isotropic quark matter takes place as crossover phase 
transition, and an increase in € can lead to a catalysis of chiral 
symmetry breaking. 


In this work, the chiral critical temperature, Tp, was de- 
termined by the peak location of the associated chiral sus- 


ceptibility Xen, which is defined as yn = amg . We stress 
that the criterion of obtaining the chiral critical temperature 
is different in different studies. Because of some shortcom- 
ings of the NJL model, such as parameter ambiguity, non- 
renormalization, and the absence of gluonic dynamics, the 
value of T, in the present work is not expected to describe 
the lattice QCD result quantitatively. However, these short- 
comings cannot affect our present qualitative results. The 
temperature dependence of chiral susceptibility Xen for dif- 
ferent € at u = 0 MeV is plotted in Fig. 1 (b). We observe 
that T, exhibits a significant é dependence. As € increases, 
T, shifts toward higher temperatures and the height of the 
peak decreases. The locations of Tę for € = —0.3, 0, 0.3 
are ~ 180 MeV, 188 MeV, and 197 MeV, respectively, which 
means a change of approximately 10% in temperature. Actu- 
ally, the in-medium meson masses also can be regarded as a 
signature of chiral phase transition. In Fig. 2(a), we display 
the variation of 7 and o meson masses with temperature for 
different € at 44 = 0 MeV. As can be observed, the 7 mass 
remains approximately constant up to a particular tempera- 
ture whereas the o mass first decreases and then increases. 
As temperature increases further, the difference between the 
m mass and o mass decreases and finally vanishes, where o 
and 7 mesons are degenerate and become nonphysical de- 
grees of freedom, which indicates the restoration of chiral 
symmetry. Further, before 7 and o meson masses emerge, 
the m mass decreases as € increase beyond T,, whereas the o 
mass first increases and then decreases with the increase in 
€. Our result slightly differs from that in Ref [40], in which 
the NJL model is used for a quark matter of finite size. The 
study shows that below the critical temperature, the 7 mass 
enhances as the system size decreases, whereas the o mass 
first remains unchanged and then increases. We also see that 
at a certain temperature, the double constituent quark mass 
(2m,) is equal to the 7 mass, so the pion meson is no longer 
a bound state but only a qq resonance and obtains a finite de- 
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Fig. 3. (Color online) Three-dimensional plot of the constituent quark mass m, with respect to temperature and quark chemical potential for 


different anisotropy parameters (€ = —0.3, 0, 0.3). 


Fig. 4. (Color online) Three-dimensional plot of chiral susceptibil- 
ity Xen for € = —0.3 in the entire u and T ranges of interest. The 
gray area means that Xen is divergent. The values remain finite due to 
numerical problems (differential quotient). The peak height at high 
H is two orders of magnitude higher than that in cases with low u 
and can be considered "divergent". 


cay width. Accordingly, the Mott transition temperature by 
the definition m (Tmot) = 2Mq(Tmot) can be obtained. The 
Mott temperatures for € = —0.3, 0, and 0.3 turn out to be ~ 
187 MeV, 196 MeV, and 206 MeV, respectively, which are 
slightly higher than the corresponding Ty. In the vicinity of 
the Mott temperature, the o meson features its minimal mass. 
In Fig. 2 (b), we illustrate the variation in the decay widths of 
both o and 7 mesons with temperature for different €. As can 
be observed, the decay width of the o meson, I's, is finite in 
the entire temperature range whereas the decay width of the 7 
meson, I+, starts after the Mott temperature. At high temper- 
ature, the merging behaviors of the decay widths for different 
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Fig. 5. (Color online) Chiral phase diagram for different anisotropy 
parameters in the Nambu—Jona—Lasinio (NJL) model. The inside 
curve is for € = —0.3, the next curve is for € = 0, and the outermost 
curve is for € = 0.3. The solid lines denote the first-order phase 
transition curves, the dashed lines denote the crossover transition 
curves, and the solid dots represent the critical endpoints (CEPs). 
We observe that the CEP is shifted towards larger values of the quark 
chemical potential but smaller values of the temperature for higher 
anisotropy parameters. 


mesons are also observed. With the increase in €, the decay 
widths of mesons are reduced. 

We continue the analysis in the finite quark chemical po- 
tential case to investigate the effect of momentum anisotropy 
on the phase boundary and CEP position. First, we display 
the temperature- and quark chemical potential-dependence of 
constituent quark mass mg for different anisotropy parame- 
ters, as shown in Fig. 3. We can observe that at a small u, 
mM, continuously decreases with increasing T, whereas mq 
has a significant discontinuity or a sharp drop along the T- 
axis at sufficiently high u, which is usually considered as 
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Fig. 6. (Color online) (a) Total cross section of the quark—quark 
scattering processes Oqq as a function of temperature at y = 0 MeV 
for different anisotropy parameters. (b) The total cross section 
of quark-antiquark scattering processes Ogg as a function of tem- 
perature at 4 = 0 MeV for different anisotropy parameters, i.e., 
€ = —0.3 (orange long dashed line), € = 0 (blue dashed line), 
and € = 0.3 (red solid line). The gray vertical lines (from left to 
right) represent the critical temperatures Te = 180 MeV, 188 MeV, 
and 197 MeV for € = —0.3, 0, 0.3, respectively. 


the appearance of a first-order phase transition. To visu- 
alize the phase diagram, we use the significant divergency 
of Xen at sufficiently high chemical potential as the crite- 
rion for a first-order phase transition, as shown in Fig. 4. 
With the decrease in u, the first-order phase transition ter- 
minates at a CEP, where the phase transition is expected to 
be of second order. As u decreases further, the maximum 
of the chiral susceptibility (<n) as the crossover criterion. 
The full chiral boundary lines in the (u-T) plane for three 
different values of € are displayed in Fig. 5. We observe 
that as € increases, the phase boundary shifts toward higher 
quark chemical potentials and higher temperatures. We ob- 
serve that the CEP appears in the low temperature and high 
chemical potential regions. Once the effect of momentum 
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Fig. 7. (Color online) Relaxation time of quarks at u = 0 MeV 


as a function of temperature for different anisotropy parameters, i.e., 
€ = —0.3 (orange long dashed line), € = O (blue dashed line), 
and € = 0.3 (red solid line). 


anisotropy (€ > 0) is taken into account, the rapid expansion 
and fast cooling of the created fireball along the beam direc- 
tion make the temperature of the anisotropic system lower 
than that of the isotropic system under the same conditions. 
By modifying the distribution function in the gap equation, 
we can study the influence of momentum anisotropy on the 
position of the CEP. It is noted that the criteria for deter- 
mining the CEP position remain the same for both isotopic 
and anisotropic systems. Therefore, as € increases, the mo- 
mentum component (temperature) in the anisotropic distri- 
bution function f'®°(4/p? + ¿(p - n)?) increases (decreases). 
Accordingly, the chemical potential in the anisotropic func- 
tion for determining the CEP position is greater than that 
in the isotropic distribution function. The CEP locations 
(Tcep, cep) in this work are presented at (298.70 MeV, 
88.2 MeV), (321.8 MeV, 82.4 MeV), (348.4 MeV, 
and 74.2 MeV) for € = —0.3, 0, 0.3, respectively. The po- 
sition of the CEP for € = 0 in this work is almost consistent 
with the existing result [118] for the same parameter set. The 
value of cep (Tcrp) from € = —0.3 to € = 0.3 increases 
(decreases) by approximately 16% (17%), indicating that the 
influence degree of momentum anisotropy on the temperature 
of the CEP is almost the same as that on the quark chemical 
potential of the CEP. This is different from the effect due to 
the finite volume in Ref. [38], which has indicated that in the 
PNJL model, the finite volume affects the CEP shift along the 
temperature stronger than along the quark chemical potential 
shift. When the system size is reduced to 2 fm, the CEP in 
the PNJL model vanishes and the whole chiral phase bound- 
ary becomes a crossover curve. Based on this result, there 
also exists a possibility that if € further increases, the CEP 
may disappear from the phase diagram. 

To better understand the qualitative behavior of the trans- 
port coefficients in the quark matter, we first discuss the 
results of the scattering cross sections and the relaxation 
time. In Fig 6, we display the total cross section of quark— 
quark scattering processes Ogg = Guu—suu + Tud>ud (plot 
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Fig. 8. | (Color online) Temperature dependence of scaled shear 
viscosity n/T’ in quark matter at vanishing chemical potential for 
different anisotropy parameters, i.e., € = —0.3 (orange long dashed 
line), € = O (blue dashed line), and € = 0.3 (red solid line). The 
thick cyan dotted line represents the result in the Ny = 3 quasiparti- 
cle model (QPM) [59], which is an effective model for the descrip- 
tion of non-perturbative QCD. The purple dash-dotted line shows 
the result obtained in the Ny = 2 NJL model by Zhuang et al. [57]. 
The brown dots show the result from hadron resonance gas (HRG) 
model [70]. The green dots correspond to the result of Rehberg et al. 
in the Ne = 3 NJL model [58] using the averaged transition rate 
method for the estimation of relaxation time. 


a) and the total cross section of quark—antiquark processes 
Og = Cunsua + Cudsud + Sua—ag (plot b) as functions of 
temperature at different anisotropy parameters for vanishing 
quark chemical potential. As can be observed, Og, and Ogg 
have similar peak structures in their temperature dependence, 
i.e., the scattering cross sections first increase, reach a peak, 
and decrease with increasing temperature afterwards. In ad- 
dition, the magnitude of Ogg is higher than that of ,,, and 
this is mainly due to the additional s-channel contribution to 
a resonance of the exchanged meson with the incoming quark 
and antiquark, which leads to a large peak in the cross sec- 
tion [74]. We further observe that the scattering cross sec- 
tions in the weakly anisotropic medium have the same behav- 
iors as those in the isotropic medium. As € increases, Ggq in- 
creases in the entire considered temperature domain, whereas 
Oqq first decreases as € increases, and then increases as € in- 
creases. With an increase in £, the maximum of the scat- 
tering cross section shifts toward higher temperatures. The 
location of the maximum for Gy, at different € is nearly in 
agreement with the respective T, as the peak positions of gg 
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Fig. 9. (Color online) Temperature dependence of scaled electri- 


cal conductivity o/T in quark matter at vanishing chemical po- 
tential for different anisotropy parameters, i.e., € = —0.3 (orange 
long dashed line), € = O (blue dashed line), and € = 0.3 (red 
solid line). The green dotted line shows the result of Marty et al. 
in the Ne = 3 NJL model [54]. The thick gray dash-dotted line 
represents the result from the pQCD-based microscopic Boltzmann 
approach to multi-parton scatterings (BAMPS) transport model [67] 
with running coupling constant. The brown stars present the result in 
the parton-hadron-string dynamics (PHSD) transport approach [66]. 
The cyan dotdashed line shows the result within the excluded vol- 
ume hadron resonance gas (EVHRG) model with the RTA [71]. The 
darkyellow dots are the lattice date obtained from Ref. [69]. The 
red open circles are the calculation for hadronic gas in the transport 
approach simulating many accelerated strongly-interacting hadrons 
(SMASH) [65] based on the Green—Kubo formalism. 


respectively locate at ~ 1.07 T7®3, 1.07 T°, 1.10 T? for 


€ = —0.3, 0, 0.3, with TÈ denoting the chiral critical tem- 
perature for a fixed €. 

The dependence of total quark relaxation time 7, on tem- 
perature for vanishing quark chemical potential at different € 
is displayed in Fig. 7. As can be observed, 7, first decreases 
sharply with increasing temperature, and after an inflection 
point (viz, the peak position of Ogg), T, changes modestly 
with temperature. Further, the increase in 7, with € is signif- 
icant at low temperature, whereas at high temperature the re- 
duction in Tg with € is imperceptible. This is the result of the 
competition between the quark number density and the total 
scattering cross section in Eq. (57). At low temperature, the € 
dependence of 7, is mainly determined by the inverse quark 
number density, whereas at high temperature it is primarily 
governed by the inverse total cross section, even though this 
effect is largely cancelled out by the inverse quark density 
effect. 
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Fig. 10. (Color online) Temperature dependence of the Seebeck co- 
efficient in quark matter at js = 100 MeV for different anisotropy 
parameters, i.e., € = —0.3 (orange broad dashed line), € = 0 (blue 
dashed line), and € = 0.8 (red solid line). The brown dotted line 
corresponds to the result for the QGP in the quasiparticle model [83] 
at {tg = 50 MeV. The cyan thick-dotted line represents the result in 
the hadron resonance gas model for ug = 0.1 GeV [81]. The mauve 
dash-dotted line and green dots, respectively, represent the results in 
the HRG model for ug = 50 MeV [80] and the Ny = 2 NJL model 
for = 100 MeV [84], where the gradient of the quark chemical po- 
tential apart from a spatial gradient in temperature is also included. 


Next, we discuss the results regarding various transport co- 
efficients. In Fig. 8, the temperature dependence of scaled 
shear viscosity n/T? in quark matter for different momen- 
tum anisotropy parameters at a vanishing chemical potential 
is displayed. We observe that with increasing temperature, 
n/T? first decreases, reaches a minimum around the critical 
temperature, and increases afterwards. The temperature posi- 
tion for the minimum n/T? is consistent with the temperature 
for the peak @,q. This dip structure of n/T? mainly depends 
on the competition between the quark distribution function f 
and the quark relaxation time 7, in the integrand of Eq. (44). 
The decreasing feature of n/T? in the low temperature do- 
main is governed by 7,, whereas in the high temperature do- 
main, the increasing behavior of ci overwhelms the decreas- 
ing behavior of 74, resulting in n/T? as an increasing function 
of temperature. In addition, we observe that with an increase 
in £, n/T? has an overall enhancement and the minimum of 
n/T? shifts to higher temperatures. This behavior of n/T? 
can be understood from the related expression in Eq. (44), 
where apart from the €-dependent relaxation time, the first 
term in the integrand of Eq. (44) has an additional € factor, 
which leads to an overall enhancement of the absolute first 
term at low T. The variation of the first term at low T is larger 
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than the counterpart of the second term, which results in a 
suppression of 7/T® for the inclusion of positive €. Mean- 
while, at high T, the qualitative and quantitative behavior of 
n/T? with £ is dominated by the second term. The location 
of the minimum for n/T? at different £ is consistent with the 
peak position of G4q. We further observe that n/T? decreases 
as € increases in the entire temperature region. In addition, 
we also compare our result for € = 0 with the results reported 
in the literature. The calculation of n/T? in the hadron reso- 
nance gas (HRG) model [70] (brown dots) using the RTA is a 
decreasing function with temperature, which is qualitatively 
similar to ours below the critical temperature. The quantita- 
tive difference between the HRG model result and ours can be 
attributed to the use of different degrees of freedom and scat- 
tering cross sections. The result of Zhuang et al [57] in the 
Ns = 2 NJL model (purple dash-dotted line) is of the same 
order of magnitude as ours, whereas at high temperature, their 
result still has a decreasing feature because an ultraviolet cut- 
off is used in all momentum integrals whether the tempera- 
ture is finite or zero. The result estimated in the quasiparticle 
(QPM) [59] is a logarithmically increasing function of tem- 
perature beyond the critical temperature, and it is quantitively 
larger than ours beyond the critical temperature owing to the 
differences in both the effective quark mass and relaxation 
time. The result of Rehberg et al [58] for the Ne = 3 NJL 
model in the temperature regime close to the critical temper- 
ature is smaller than ours, and the obvious dip structure is not 
observed because the momentum cutoff is also used at finite 
temperature. 


In Fig. 9, we plot the thermal behavior of scaled electrical 
conductivity o./T at y = 0 MeV for different €. Similar to 
the temperature dependence of 7/T°, oe|/T also exhibits a 
dip structure in the entire temperature region of interest. We 
also present the comparison with other previous results. The 
result obtained from the PHSD approach [66] (brown stars), 
where the plasma evolution is solved by a Kadanoff—Baym 
type equation, also has a valley structure, even though the lo- 
cation of the minimum is different from ours. We also observe 
that in the temperature region dominated by the hadronic 
phase, the thermal behavior of o./T using the microscopic 
simulation code SMASH [65] (pink open circles) is similar to 
ours. Furthermore, our result is much larger than the lattice 
QCD data (dark yellow dots) taken from Ref. [69] owing to 
the uncertainty in the parameter set and absence of gluonic 
dynamics. The result within the excluded volume hadron res- 
onance gas (EVHRG) model [71] (cyan dash-dotted line) and 
the result obtained from the partonic cascade BAMPS [67] 
(gray thick dash-dotted line) are qualitatively and quantitively 
similar to our calculations below the critical temperature and 
beyond the critical temperature, respectively. Our result is 
similar to that of Marty et al. obtained within the Nr = 3 
NJL model [54] (green dotted line), with the numerical dis- 
crepancy mainly coming from the differences in the values of 
the model parameter set and scattering cross sections. At low 
T, the absolute values of both the first and second terms in 
Eq. (54) increase as € increases. However, the variation of 
the first term is larger than that of the second term, which re- 
sults in an enhancement of ca /T. At high T, the decreasing 


feature of relaxation time with € can weaken the increasing 
behavior of ce /T with £, and the values of a. /T for different 
€ gradually approach and eventually overlap. Our qualitative 
result of o./T is different from the result in Ref. [73], where 
the o./T of the QGP is a monotonic decreasing function of £. 
This occurs because the effect of momentum anisotropy is not 
incorporated in the calculation of the relaxation time and the 
effective mass of quasiparticles, as the € dependence of oa /T 
is only determined by the anisotropic distribution function. 
We also observe that with the increase in €, the minimum of 
oa/T shifts to higher temperatures, which is similar to 7/T®. 
However, the height of the minimum increases, which is op- 
posite to 7/T®. 


Finally, we study the Seebeck coefficient S in quark- 
antiquark matter. Owing to the sensitivity of S to the charge 
of particle species, at a vanishing chemical potential, quark 
number density n, is equal to antiquark number density ng, 
and the contribution of quarks to S' is exactly cancelled by 
that of antiquarks. Thus, a finite quark chemical potential 
is required to obtain a non-zero thermoelectric current in the 
medium. In Fig. 10, we plot the variation of S with respect to 
temperature for different € at u = 100 MeV. The comparison 
with other previous calculations, which were all performed 
in the kinetic theory under the RTA, is also presented. We 
remind the reader that at finite 4, ng is larger than ng, and 
the contribution of quarks to total S in magnitude is always 
prominent. As shown in Fig. 10, the sign of S' in our investi- 
gation is positive, which indicates that the dominant carriers 
converting the heat gradient to the electric field are positively 
charged quarks, i.e., up quarks. Actually, the positive or neg- 
ative sign of S' is mainly determined by the factor (Eq — Hq) 
in the integrand of Eq. (55). In Ref. [83], the Seebeck coeffi- 
cient studied in the QPM (brown dotted line) at 4 = 50 MeV 
also exhibits a decreasing feature with increasing tempera- 
ture. The result of Abhishek et al [84] at y = 100 MeV 
in the Ne = 2 NJL model (the green dots) is very different 
from ours. In Ref. [84], S is negative and its absolute value is 
an increasing function with temperature. The reasons behind 
this quantitative and qualitative discrepancy are twofold: (1) 
the relaxation time in Ref. [84] was estimated using the av- 
eraged transition rate w;;, whereas our relaxation time was 
obtained from the thermally averaged cross section of elastic 
scattering (a detailed comparison of the two methods can be 
found in Ref. [74]); (2) in Ref. [84], the spatial gradient of 
chemical potential was also included apart from the temper- 
ature gradient and accordingly, the sign of S was primarily 
determined by the factor (E, — w/nq) with w denoting the 
enthalpy density in the associated formalism. Given that the 
single-particle energy Æq remains smaller than (w/n,), S in 
Ref. [84] is negative. We also observed that with increasing 
temperature, S sharply decreased below T., whereas the de- 
creasing feature of S was inconspicuous above T}. Further, 
the value of S at low T was much larger than that at high T. 
This is also different from the result in Ref. [84], where the 
absolute value of S in quark matter increased with increasing 
temperature because of the increasing behaviors of both the 
factor | — w/n,| and the equilibrium distribution function. In 
addition, the Seebeck coefficient in the HRG model [80, 81] 
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is also positive (negative) without (with) the spatial gradient 
of chemical potential (cyan thick dotted line and mauve dash- 
dotted line). Nevertheless, the absolute value of S in hadronic 
matter is still an increasing function of temperature regardless 
of the spatial gradient of u. We also observed that as £ in- 
creases, S has a quantitative enhancement, which is primarily 
due to a significant rise in the thermoelectric conductivity a, 
even though 1/c« has a cancellation effect on the increase in 
S. At sufficiently high temperature, the rise in 1/o,, can al- 
most compensate the reduction in a, and as a result, S varies 
insignificantly with the € of interest, compared to the value of 
S itself. 


VII. SUMMARY 


We phenomenologically investigated the impact of weak 
momentum-space anisotropy on the chiral phase structure, 
mesonic properties, and transport properties of quark matter 
in the 2-flavor NJL model. The momentum anisotropy, which 
is induced by the initial preferential expansion of the created 
fireball in heavy-ion collisions along the beam direction, can 
be incorporated in the calculation through the parameteriza- 
tion of the anisotropic distribution function. Our result has 
shown that the chiral phase transition is a smooth crossover 
for vanishing quark chemical potential, independent of the 
anisotropy parameter €, and an increase in € can even hinder 
the restoration of the chiral symmetry. We found that the CEP 
is highly sensitive to the change in €. With the increase in £é, 
the CEP shifts to higher u and smaller T, and the momentum 
anisotropy affects the CEP temperature to almost the same 
degree as it affects the CEP chemical potential. Before the 
merge of m and o meson masses, the € dependence of the 7 
meson mass is opposite to that of the a meson mass. 

We also studied the thermal behavior of various transport 
coefficients, such as the scaled shear viscosity n/T, scaled 
electrical conductivity o./T, and Seebeck coefficient S at 
different €. The associated £-dependent expressions are de- 
rived by solving the relativistic Boltzmann—Vlasov transport 
equation in the relaxation time approximation, and the mo- 
mentum anisotropy effect is also embedded in the estimate 
of relaxation time. We found that 7/T? and o./T have a 
dip structure around the critical temperature. Within the con- 
sideration of momentum anisotropy, n/T? decreases as € in- 
creases and the minimum shifts to higher temperatures. With 
the increase in £, o./T significantly increases at low tem- 
perature, whereas its sensitivity to € at high temperature is 
significantly reduced, which is different from the behavior of 
n/T? with respect to €. We also found that the sign of S$ at 
u = 100 MeV was positive, indicating that the dominant car- 
riers for converting the thermal gradient to the electric field 
are up quarks. With increasing temperature, S first decreases 
sharply and then almost flattens out. At low temperature, S 
significantly increases with the increase in €, whereas at high 
temperature the rise is marginal compared to the value of S 
itself. 

We note that it is of considerable interest to include the 
Polyakov-loop potential in the present model to study both 


chiral and confining dynamics in a weakly anisotropic quark 
matter. A more general ellipsoidal momentum anisotropy 
characterized by two independent anisotropy parameters is 
then needed to gain a deeper understanding of the QGP prop- 
erties. In the present work, no proper time dependence was 
given to the anisotropy parameter. However, in a realistic 
case, € varies with the proper time starting from the initial 
proper time up to a time when the system becomes isotropic. 
Thus, a proper time dependent anisotropy parameter [119] 
needs to be introduced to better explore the effect of time- 
dependent momentum anisotropy on chiral phase transition. 
For the strongly longitudinal expanding QCD matter, the in- 
vestigation of chiral phase transition needs to be performed by 
numerically solving both the Vlasov equation and gap equa- 
tion concurrently and continuously. In this case, the phase di- 
agram of a strongly expanding system is a map in the space— 
time plane rather than in the T — u plane. In addition, the 
investigation of the thermoelectric coefficients, especially the 
magneto-Seebeck coefficient and Nernst coefficient in mag- 
netized quark matter, based on the PNJL model would be an 
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attractive direction, and we plan to work on it in the near fu- 
ture. 


APPENDIX 


In the Nf = 2 NJL model, there are 12 different elastic 
scattering processes: 


uli > uti, ud > ud, uti > dd, 
uu > uu, ud > ud, uti > UU, 
ud — td, dd > dd, dd > uŭ, 

du — di, dd > dd, dd > dd. (62) 


The explicit expressions of the squared matrix elements for 
the uŭ > ut, ud — ud, and ud > ud processes via ex- 
change of scalar and/or pseudoscalar mesons to 1/N, order 
are given as 


2 1 g < 
|Mua—+ual?(s,t) = s?|D7|? + t|D7|? + (s — 4m?)?|D?|? + (t — 4m’)? | DF? + wie jan: Dt + s(4m° —t)D™ D? 


+t(4m? — s)D™D? + (st + 4m?(s +t) — 16m") Dz De | 


|Mud-+ual’(8,t) = 48| DF? + t?|DF/? + (¢ 


|Mud>ual (t, u) = 4u?|D™|? + | DT /? + (t 


The meson propagators in the above processes are €- 
dependent. Based on the above formulae of three scattering 
processes, the squared matrix element for the remaining scat- 
tering processes can be obtained through charge conjugation 
and crossing symmetry [57, 120]. 
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